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1, INTRODUCTION 

Nowadays, liquid slosh inside a cargo always happens in many situations. For example, ships with 
liquid container carriers are at high risk of generating sloshing load during operation [1]. In the metal 
industries, high oscillation can spill molten metal that is dangerous to the operator [2]. Meanwhile, sloshing 
of fuel and other liquids in moving vehicles may cause instability and undesired dynamics [3]. Hence, it is 
necessary to completely study the behavior of this residual slosh induced by the container motion. One may 
study the behavior of liquid slosh through developing the exact mathematical model of liquid slosh. So far, 
many researchers focus on the first principle approach to model the slosh behavior, while there are few 
literatures to discuss it from the perspective of nonlinear system identification approach. 

On the other hand, block oriented nonlinear system identification has became a popular techniques 
to model a complex plants. The block oriented nonlinear model can be classified into three categories, which 
are Hammerstein model, Wiener model and Hammerstein Wiener model. In particular, Hammerstein model 
is a model that consists of a nonlinear function followed by linear dynamic sub-plant, while Wiener model 
consists of a linear dynamic sub-plant followed by nonlinear function, and finally, Hammerstein-Wiener 
model contains a linear dynamic sub-plant inserted between two or more nonlinear functions in series. 
Among these three block oriented models, Hammerstein model is famous due to its simple model structure 
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and it has been widely used for nonlinear system identification. Specifically, the Hammerstein model has 
been applied to model a real plant such as Solid Oxide fuel cell [4], bidirectional DC motor [5], oxygen 
uptake estimation [6], stretch reflex dynamics [7], turntable servo system [8], pneumatic muscle actuators [9], 
amplified piezoelectric actuators [10] and multi-axis piezoelectric micro positioning stages [11]. On the other 
hand, there are many tools that have been utilized to identify the Hammerstein model. There are the iterative 
method [12-14], the subspace method [15-17], the least square method [18], the blind approach [19] and 
the parametric instrumental variables method [20]. Moreover, many also consider the optimization tools for 
Hammerstein model, such as Bacterial Foraging algorithm [21], Cuckoo search algorithm [22], Particle 
Swarm optimization [23], and Genetic algorithm [24]. Based on the above literature, several limitations are 
ineluctable in their works, which are: 
a. Most of the Hammerstein models used in their study are based on discrete-time model, while many real 
plants can be easily represented in continuous-time model. 
b. Almost all the methods assume a known structure of nonlinear function, which consists of several 
basis functions. 

Though, our proposed work can solve a more general class of continuous-time Hammerstein model by 
assuming an unknown structure of nonlinear function. In particular, a piece-wise affine function is adopted 
with so many basis functions. Due to the introduction of the piece-wise affine function, a high dimensional 
design parameter tuning is considered in this study, which make the identification problem more complex. 
On the other hand, Grey Wolf Optimizer (GWO) [25] has become a top notch optimization algorithm which 
has solved various types of engineering problems [25-27]. To the best of our knowledge, there are still few 
works to discuss on the GWO for identification of Hammerstein model. Thence, it motivates us to see 
the effectiveness of the GWO in modeling the liquid slosh plant from the real experimental data. 

This paper presents the identification of liquid slosh plant using the Hammerstein model based on 
GWO method. A remote car carrying a container of liquid is considered as the liquid slosh experimental rig. 
The GWO method is used to tune both coefficients in the nonlinear function and transfer function of 
the Hammerstein model such that the error between the identified output and the real experimental output is 
minimized. The effectiveness of the proposed framework is assessed in terms of the convergence curve 
response, output response, and the stability of the identified model through the bode plot and pole zero map. 


2. RESEARCH METHOD 

This section explains on how the identification of liquid slosh based on Hammerstein model and 
GWO based method is performed. Firstly, the experimental setup of liquid slosh plant is explained. Here, it is 
shown on how the input and output data are generated. Secondly, we show the procedure to identify 
the model of liquid slosh plant using the Hammerstein based Grey Wolf Optimizer. 


2.1. Liquid slosh experimental rig 

In this study, a mobile liquid slosh plant is considered to replicate real situation of a moving 
container carrying liquid, as shown in Figure |. In particular, a remote control car is used to carry a small 
tank filled with liquid. The tank is also equipped with four plastic wheels so that it can move smoothly as 
shown in Figure 1(a). Moreover, three accelerometer sensors (ADXL335) that are floated on the surface of 
liquid are used to measure liquid oscillation as shown in Figure 1(b). For simplicity of our study, the liquid 
slosh data from only one of the sensor is recorded and only z-axis output data is considered. Figure 2 shows 
a general schematic diagram of liquid slosh experimental rig. 





(a) (b) 
Figure 1. Liquid slosh experimental rig, (a) side view, (b) plan view 


Grey wolf optimizer for identification of liquid slosh behavior... (Mohd Ashraf Ahmad) 


544 0 ISSN: 2302-9285 





Monitor —— Liquid Slosh plant 


Data Acquisition 


Figure 2. Schematic diagram of liquid slosh experimental rig 


In particular, an Arduino UNO 1s used as a data acquisition platform to process the input and output 
data. Here, we generate a voltage from the Arduino UNO to the remote car and concurrently the Arduino 
UNO also will acquire the slosh data from the accelerometer. Both the input and output data can be 
monitored and analyzed from the personal computer using the LabView software. In order to identify 
the model of liquid slosh, the remote car is required to move to a certain distance and suddenly stop to 
generate a liquid oscillation or slosh inside the tank. Thence, we apply the input voltage as shown in Figure 3 
to move the remote car. Concurrently, the liquid slosh data is recorded as shown in Figure 4. These two data 
are then used to develop the Hammerstein model based GWO, which is discussed in the next section. 
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Figure 3. Input voltage applied to the remote car Figure 4. Output slosh from the accelerometer 


2.2. Identification of liquid slosh using Hammerstein based grey wolf optimizer 

In this section, the proposed Grey Wolf Optimizer (GWO) for identification of liquid slosh plant in 
Section 2.1 based on Hammerstein model is presented. Firstly, a problem formulation to identify the liquid 
slosh plant is explained. Then, it is shown on how to apply the GWO method to identify the liquid slosh 
based on Hammerstein model. Figure 5 shows a complete block diagram to identify the liquid slosh model in 
Section 2.1. The proposed Hammerstein model consists of nonlinear function h(u) followed by the transfer 
function G(s). The nonlinear function is a piece-wise affine function given by 





Figure 5. Block diagram of Hammerstein model based GWO 
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Co + m,(u — do) if dgsu<d, 
= i < 
h(u) = Cy + M>(u d,) | if a, su< d» (1) 
Co-1t+mo(u-—d,_-1) if dg .~<u<d, 
and the transfer function G(s) is given by 
B(s s™ +b, 4s" 1 +--+ 
G(s) = ( ) _ m-1 0 (2) 





~— A(s) — AmS™ + An_1S™ 1 + + Ag 


In (1), the symbol m; = (c; — c;_1)/(d; — dj_,)(@i = 1,2, ...0) are the segment slope with connecting input 
and output points as d;(i=0,1,...0)and c;(i =0,1,..0), respectively. For simplicity of notation, 
let d=[do, di, ..., ds]! and c=[co, C1, ..., Co]'. The input of the real liquid slosh plant and the identified model is 
defined by u(t), while the output of the real liquid slosh plant and the identified model are denoted by y(t) 
and 7(t), respectively. Thence, the expression of the identified output can be written as 


H(t) = G(s)h(u(t)) (3) 


Moreover, several assumptions are adopted in this work, which are: 
a. The order of the polynomial A(s) and B(s) are assumed to be known 
b. The nonlinear function A(u(f)) is one-to-one map to the input u(t) and the values of d;(i = 1,2, ...0) are 
pre-determined according to the response of input u(f). 
Next, let t; be a sampling time for the real experimental input and output data (u(t), y(t)) 
(t=0, ts, 2ts, ..., Nts). Then, in order to accurately identify the liquid slosh model, the following objective 
function in (4) is adopted in this study: 


EG,h) = > (vats) - IGts))” (4) 
n=0 


Finally, our problem formulation can be described as follows. 

Problem 1. Based on the given real experimental data (u(t), y(t)) in Figure 1, find the nonlinear 
function h(u) and the transfer function G(s) such that the objective function in (4) is minimized. 

Furthermore, it is shown on how to apply the GWO in solving Problem 1. For simplicity, let the 
design parameter of Problem 1 is defined asx=[bp by .. Dy-1 GQ UG += An Co = Col’, 
where the elements of the design parameter are the coefficients of both the nonlinear function and the transfer 
function of the continuous-time Hammerstein model. In GWO framework, let x;(i = 1,2,...,M) be the 
design parameter of each agent i for M total number of agents. Then, consider x;;(j = 1,2,...,D) be the j-th 
element of the vector x;(i = 1,2,...,M), where D is the size of the design parameter. Thence, by adopting 
objective function in (4), a minimization problem is expressed as 


OF xiao EQ) (9) 

for iterations k=1, 2, ..., until maximum iteration kmax. Finally, the procedure of the GWO in solving 

Problem 1 is shown below: 

- Step 1: Determine the total number of agents M and the maximum iteration kmax. Set k=O and initialize 
the design parameter x;(0)(i = 1,2, ...,M) according to the upper bound Xup and lower bound Xow values 
of the design parameter. 

- Step 2: Calculate the objective function in (4) for each search agent i. 

- Step 3: Update the values of the best design parameter x*, the second best design parameter x’ and 
the third best design parameter x’ based on the generated objective function in Step 2. 

- Step 4: For each agent, update the design parameter using the following equation: 


xi +x%7+%x3 
x(k + 1) 7 (6) 


where 


xt =x%—P-|Q-x% —x,(k)|,x? =x? —P-|Q- x8 —x,(k)|, 


x3 = x° —P-|Q-x° —x({k)I ) 
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for i=1, 2, ..., M. In (6), the vectors P and Q are expressed as follows: 
P=2a-:rn-—a (3) 
Q=2:'7r (9) 


where r; and 72 are the random vectors, where each element is generated independently from 0 to 1, and each 
element of a is linearly decreased from 2 to 0 over the course of iteration using the following equation: 





a=2(1- : (10) 


max 


Note that the vectors of P and Q in (8) and (9), which consists of random vectors r; and r2, are generated 

independently between x!, x? and x°. 

- Step 5: After the maximum iteration is achieved, record the best design parameter x* and obtained 
the continuous-time Hammerstein model in Figure 1. Otherwise, repeat Step 2. 


3. RESULTSAND ANALYSIS 

In this section,the effectiveness of the GWO based method for identifying the liquid slosh system 
using continuous-time Hammerstein model is demonstrated. In particular, the convergence curve response 
of the objective function in (4), the bode plot and pole-zero mapping of linear function and the plot 
of nonlinear function, will be presented and analyzed in this study. 

Based on the experimental setup in Section 2.1, the input response u(t) as shown in Figure 3 
is applied to the liquid slosh plant, and the output response y(f) is recorded as shown in Figure 4. 
Here, the input and output data are sampled at 4,=0.02 for N=450. In this study, the structure of G(s) is 
selected as follows: 


B(s) _ s? + bys* + bys + bo 


G(s) = —— —————— ————————————————— 
(s) A(s) a,s*++a 3s? + a,s*+a,5+ a) 


(11) 





After performing several preliminary testing on the given data (u(t), y(t)). Meanwhile, the input points 
for piece-wise affine function of A(u(t)) are given by d=[0, 0.2, 0.4, 0.6, 0.8, 1, 2, 3, 4, 5]'. The design 
parameter x € R?® with its corresponding transfer function and nonlinear function is shown in Table 1. 
Next, the GWO algorithm is applied to tune the design parameter with initial values of design parameter are 
randomly selected between the upper bound Xyp and lower bound Xow as shown in Table |. Here, we choose 
the number of agents M=40 with maximum iterations Kmax=5000. 


Table 1. Design parameter of liquid slosh plant 


x Coefficients Xiow Xup x 

x1 b> —5 35 —4.8017 
X2 by —5 35 0.2657 
X3 bo —5 35 —0.0243 
XA a4 —5 35 —0.1119 
x5 a3 —2200 —1 —4,3214 
X5 a2 —2200 —| —31.6384 
x7 a) —2200 —| —901.2917 
Xe ao —2200 —1 —1280.4483 
Xo Co —5 5 —1.9997 
X10 C\ —5 5 0.3249 
xu C2 —5 5 1.5968 
X12 C3 —5 5 —0.0757 
X13 C4 —5 ) —0.0099 
X14 C5 —5 5 —0.0048 
X15 C6 —5 5 —0.0251 
X16 C7 —5 5 0.0584 
X17 Cg —5 5 0.0216 
X18 C9 —5 5 —2.2603 
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Figure 6 shows the response of the objective function convergence after 5000 iterations to produce 
the best design parameter as shown in the final column of Table 1. It shows that the GWO based method 
is able to minimize the objective function in (4) and produce a quite close output response y(t) as compared 
to the real output y (t), which can be clearly seen in Figure 7. Note that the identified output response tends to 
yield high oscillation when input is injected to the system and it start to attenuate when the input is zero, 
which is quite similar to the response of real experimental output. 
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Figure 6. Convergence curve response Figure 7. Response of the identified output (t) 
and real output y(t) 


In the real experimental setup, we can say that the liquid slosh system is stable since the liquid slosh 
output is reduced gradually as t — oo. In order to validate our model regarding the stability, we use the bode 
diagram and the pole-zero map of the identified transfer function G(s) as shown in Figures 8 and 9, 
respectively. From the bode diagram, it can be observed that a gain margin of 2.56 dB at 17.4 rad/s and phase 
margin of 4.18 degree at 16.5 rad/s are produced, which indicates that the linear system is stable. This is also 
supported by the pole-zero map where all the poles are located at the left hand side of y-axis. In particular, 
the obtained values of poles are -36.6623, -0.2344+j14.5053 and -1.4828, while the obtained values 
of zeros are 4.7468, 0.0274+70.0661. On the other hand, we also can observe the feature of nonlinear function 
by plotting the obtained piece-wise function as depicted in Figure 10. Note that our nonlinear function is not 
restricted to any form of nonlinear function (i.e., quadratic), which is more generalized and provide more 
flexibility of searching a justifiable function. 


Bode Diagram 
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Figure 10. Resultant of piece-wise affine function h(u) 


4. CONCLUSION 

In this paper, an identification of liquid slosh plant using continuous-time Hammerstein model based 
on Grey Wolf Optimizer (GWO) has been presented. The results demonstrated that the proposed generic 
Hammerstein model based on GWO has a good potential in identifying the real liquid slosh behavior. 
In particular, it is shown that the proposed method is able to produce a quite close identified output with real 
liquid slosh output. Moreover, the resultant linear model has been proved to be stable based on the bode plot 
and pole-zero map. It is also shown that the used of piecewise-affine function gives more flexibility 
for the GWO to search more generic nonlinear function. In the future, our work can be extended to various 
types of nonlinear function such as continuous-time Wiener and Hammerstein- Wiener. 
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